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Abstract 



We compute the finite temperature pole mass shifts of the octet and de- 
cuplet baryons using heavy baryon chiral perturbation theory and the 1/Nc 
expansion, where Nc is the number of QCD colors. We consider the temper- 
atures of the order of the pion mass m^r , and expand truncate the chiral and 
l/Nc expansions assuming that ~ ^/Nc- There are three scales in the 
problem: the temperature T, the pion mass m^r, and the octet-decuplet mass 
difference. Therefore, the result is not simply a power series in T. 

We find that the nucleon and A temperature mass shifts are opposite in 
sign, and that their mass difference changes by 20% in the temperature range 
90 MeV < T < 130 MeV, that is the range where the freeze out in relativistic 
heavy ion collisions is expected to occur. 

We argue that our results are insensitive to the neglect of 1/Nc- supressed 
effects; the main purpose of the 1/Nc expansion in this work is to justify our 
treatment of the decuplet states. 
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I. INTRODUCTION 



In this paper we compute the temperature-dependence of the pole mass of the J = 2 
and </ = f baryons using effective lagrangians and large methods. There are two main 
motivations for this work. The first comes from relativistic heavy ion collisions. These are 
very complex phenomena, and the extraction of physically interesting quantities from the 
experimentally accessible information necessarily involves a good deal of modelling. Most 
likely, a number of measurable quantities will have to be correlated to constrain the ones 
of physical interest. Clearly, an understanding of the parameters of these models within 
controlled approximations will simplify this process. The dependence of particle proper- 
ties on the temperature is one such parameter. The second motivation is more theoretical. 
The finite temperature properties of hadrons provide clues to the mechanism of confine- 
ment/chiral symmetry breaking which is expected to occur at T ~ 200 MeV. Below this 
critical temperature, we can use effective lagrangian methods to gain some information about 
the finite-temperature dependence of QCD quantities. 

Chiral perturbation theory (xPT) is the effective low-energy theory of hadrons and gives 
rise to a systematic expansion in powers of the light current quark masses mu,d,s/-^ and 
kinematic scales E/A, where A ~ 1 GeV is the chiral expansion scale. All the information 
contained in QCD is encapsulated in a few phenomenological constants, that have to be 
determined by experiment. An essential result is that at a given order in the expansion one 
has to compute only diagrams with a finite number of loops depending on a finite number 
of couplings IQ. Tree diagrams correspond to current algebra results, of which xPT can 
be viewed as a systematic extension to all orders. Formulated originally for the Goldstone 
bosons TT, K, and rj, it was later extended to include the low-lying baryons 0. A non- 
trivial aspect of this extension is the presence of a new scale, the baryon mass, that does not 
vanish in the chiral limit. This appears to spoil the scaling arguments used to get the relation 
between number of loops and powers of quark mass (or momentum). An elegant solution 
to this problem is given by heavy baryon chiral perturbation theory (HBxPT) 0, which 
treats the baryon as a static particle. Corrections suppressed by powers of the baryon mass 
appear as higher-order vertices in the effective lagrangian. Since the baryon mass appears 
only in vertices, the chiral counting manifestly goes through and the relation between the 
number of loops and powers of quark mass is recovered. 

One difficulty with HBxPT is the large number of phenomenological constants appearing 
in the effective lagrangian. Another one is the presence of A states close to the nucleons 
(Ma— Mtv — 300 MeV). There are different prescriptions for dealing with the decuplet states: 
some authors include them as explicit degrees of freedom (treating the the decuplet-octet 
mass difference either perturbatively or to all orders), while others keep only the nucleon 
fields and consider only processes at energies low enough to avoid exciting the decuplet 
states. The 1/Nc expansion (where Nc is the number of QCD colors) alleviates both of 
these problems [Q. It relates the value of some phenomenological constants, allowing one 
to go further in the chiral expansion. Also, in the 1/Nc expansion there is a whole tower 
of degenerate baryon states with spins |,|,... with mass differences between low- lying 
spin multiplets suppressed by l/Nc. The l/Nc expansion therefore forces us to include the 
decuplet as explicit degrees of freedom, and (after deciding how to book l/Nc corrections 
compared to chiral corrections) the treatment of the decuplet-octet splitting is determined. 
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Another appealing aspect of the 1/Nc expansion is that the intuitive and phenomenologically 
successful! static quark model arises as its leading term HBxPT together with the l/N,, 
expansion has been used to compute a number of quantities sucessfully . In particular, a 
discussion of the zero temperature baryon mass differences at next-to-leading order will be 
reported elsewhere 

At temperatures T well below the critical temperature ~ 200 MeV, we expect a 
description of a QCD plasma in terms of hadronic degrees of freedom to be justified. Since 
most particles will have kinetic energies of order T or smaller, their interactions will be soft 
and we expect that xPT will give an accurate description of the dynamics. Particles with 
high energy or large masses are exponentially suppressed by Boltzmann factors and can 
be safely ignored. For example, Ref. [H studied the thermodynamics of a gas of pions to 
compute the temperature dependence of the chiral condensate. 

The properties of baryons in a thermal plasma are affected by the background of thermal 
pions. In this paper we compute the leading order temperature dependence of the pole mass 
of octet and decuplet baryons using HBxPT and the l/N^ expansion. 

This paper is organized as follows. In section II we explain the power counting that 
selects the class of diagrams contributing to the leading order temperature baryon pole 
mass dependence. In Section III the results are presented and discussed. The chiral limit 
is briefly considered in Section IV where a few comparisons with previous results are made. 
Finally, in Section V the main results are summarised and comments on the range of validity 
of our approach are presented. 



II. EFFECTIVE LAGRANGIAN AND POWER COUNTING 



The calculation will involve the real time formalism (RTF) for finite temperature field 
theory and the formalism presented in for large heavy baryon chiral perturbation 
theory. Here we confine ourselves to defining the notation used, leaving the details to the 
references. 

The complications caused by the presence of a matrix structure in the propagators in 
the RTF [|I^ is offset by the fact that it is easy to sepatate the temperature-dependent 
contributions from the zero temperature part. The propagator for the (static) baryon with 
velocity (taken to be the same as the plasma velocity) is 
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where n{x) {np{x)) is the bosonic (fermionic) distribution function 
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The indices 1, 2 refer to the the matrix structure of the RTF propagator needed to describe 
the expectation value of operators in a mixed state. Note the presence of the large baryon 
mass M in the distribution function np. They appear because, in HBxPT, the momentum 
in the propagator (0) is the residual momentum k = p — Mv where p is the physical momen- 
tum 10. Thus, in the range of temperatures considered here, the temperature-dependent part 
of the baryon propagator is suppressed by a tiny Boltzmann factor and can be safely dropped. 
To a lesser extent, the same is true for K and rj mesons: since n{mT^/T)/n{mK/T) ~ 10 
at T = 200 MeV, we will count thermal kaon and eta loops as addionally suppressed in 
the chiral counting. In the present leading order calculation this means they will also be 
dropped. 

In the large- iVc world, the low lying baryons form an infinite tower of flavor multiplets 
with spin i, |, . . .. We follow the formalism and notation of Refs. P,p!2[| where they are all 
treated simultaneously by defining a field acting on a spin-flavor Fock space 
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where oi, . . . , oat = 1, • • • , Np are flavor indices and ai, . . . , =t) i are spin indices. The 
operators d^ and d are (bosonic) creation and annihilation operators. The effective la- 
grangian then has the standard form 



L = {B\{tv ■ D} + fi{a,}{a^} - {M} + 9{A,a^} 
where the braces around a generic operator stand for 

{w} = w^^pdi^d'p. 

M. is built from the quark mass matrix M 
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D\B)^{d-i{V})\B), (8) 
V,^\ {iU^ + i^U) , (9) 
A,^'- {^d,e - eO,^) , (10) 

where ^(x) parametrizes the Goldstone bosons 

^(x) = e*n(x)//^ 

and / ~ 93 MeV is the pion decay constant. 

The main result obtained in Ref. is that the coefficient of an n-body operator in the 
effective lagrangian above has a coefficient which is < 1/N^~^ for large Nc. This comes 
about because, at the QCD level, an exchange of n — 1 gluons is necessary to produce a term 
in the effective lagrangian involving n quarks. This is true even when the number of light 
quark flavors Np is taken to be of order Nc, since this rule does not rely in the suppresion of 



quark loops |T^. For Np > 3, the size of the flavor representations of the low lying baryons 
with given J ~ 1 grow with and it is not obvious which of those states should be identified 
with the states of the real = 3 world. Our procedure will be to compute the mass shifts 
keeping all terms that are of the order desired in any state in the large flavor multiplet. 
Only then those terms are evaluated with N = Np = 3. 

Because we take Np ~ N^, a decision has to be made about how to extrapolate the quark 
mass matrix to arbitrary number of flavors. We choose the most general form 
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and keep all terms that are of the desired order in any extrapolation in which Nu + Nd + Ng = 
Np. 

Now we can establish the rules for power counting. We take to be of the same 
order as a strange quark mass m^. The light quark masses mu^d are counted as ~ l/N"^. 
As mentioned above, kaon and eta thermal loops are suppressed by large Boltzman factors 
compared to pion thermal loops, and are therefore negligible in this expansion. 

This suppression however, does not become exact in the chiral or large A^ limit. To 
account for this numerical fact we take the practical attitude of counting those loops as 
additionally suppressed by a factor ~ 1/A^^. Summarizing, 

— ~ m-s ~ e, m = -{rriu + md) ^ K,r] thermal loop ~ e^. (13) 

Let us now look at the graphs giving the leading temperature dependent contributions 
to the baryon mass. First, consider the graph in Fig. 1. The tick on the meson propagator 
denotes the temperature dependent part of the propagator. This graph contains two one- 
body vertices, and two powers of / in the denominator. Since / ~ \/T^c and the matrix 
element of any n— body operator is at most ~ A"", using dimensional analysis we can see that 
this graph is ~ Nmll{mn/T), where I{x) is a some dimensionless function. For T ~ we 
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have I{mT,/T) ~ 1, so this graph is booked as Ne"^ . An exphcit calculation shows, however, 
that I{mT^/T) vanishes. 

The graphs in Fig. 2, with an arbitrary number r of mass insertions of the form /x{cr*}{cr*} 
are of the order N^~^m^~^ ~ A^e^ (since fi ~ 1/N). The sum of these graphs therefore gives 
the leading contribution to the temperature mass shift. 

Note that naively, the diagram in Fig. 2 with one insertion of the operator {Ai} is 
of order N^rrism'^. This contradicts the fact that baryon masses are proportional to A^. 
When wave function renormalization terms are included however, there is a cancellation 
that decreases this result by two powers of A^, making this contribution not only consistent 
with large A"c expectations, but actually negligible at leading order in our expansion. We 
assume that similar cancellations (necessary for the consistency of large Nc QCD) occur so 
that no contribution grows faster than ~ Ac. This has been checked in various non trivial 
cases . 

Using similar arguments, one can check that other graphs including l/A^ suppressed or 
higher derivative vertices, as well as light quark mass insertions and higher loop graphs are 
also of higher order. Thus, the leading contribution for the temperature dependent part of 
the pole mass is given by the simple one loop graph with the baryon propagator resummed to 
include the (zero temperature) octet-decuplet mass splitting. Since this graph is to be eval- 
uated at the tree level mass shell (including the effects of the term yu{cr*}{cr*}), the external 
momentum will cancel against the mass in the propagator if both the internal and external 
particles, have the same spin. When this happens the graph with the resummed propagator 
is exactly equal to the graph in Fig. 1, which vanishes. Thus, the only intermediate states 
that contribute to the octet temperature mass shifts are decuplet intermediate states, and 
vice-versa. 



III. RESULTS 



The calculation of the graph in Fig. 2 with the resummed propagator is straighforward. 
Keeping only the temperature dependent part we have 
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where T"^ are SUpi^'i) generators and P is a projection operator 
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Here is the projection operator for octet masses, which projects onto decuplet interme- 
diate states; vice versa for the decuplet masses. The upper(lower) sign in ( 0) applies to 
the octet (decuplet) and X is the operator 
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-^"8 = 15 - y {5} + liSa'Ha^} + -{S}{S}, (16) 
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is the strangeness matrix in flavor space. The expectation value of X in different states is 
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The zero value for X^ reflects the fact that the Q interactions with the background pions 
involve the annihilation of two quarks into gluons that then turn into the quark making up 
the pion, which is a 1/A^ suppressed process. 

The integral in ( |1^) is to be understood in the principal value sense. The coupling 
constant g is related to the usual octet axial couplings D and F, as well as the axial couplings 
to the decuplet states. In the large- Ai'c and SUp{3) limits we can identify as g = 3{D+F)/5 ~ 
0.8, which is also the value obtained by performing a best fit at leading order. In Fig. 3 we 
show the dependence of the real part of the nucleon pole mass as a function of temperature. 
(The temperature dependence dependence for other octet baryons is a multiple of the same 
graph.) A striking feature is that it is very small until T ~ 80 MeV, and then varies rapidly. 

The shift in the mass of the A is opposite to the nucleon, so that the mass splitting 
between then decreases rapidly in the range 90 Mev < T < 130 Mev. This fact may be used 
to predict relative yields of different baryon species in central rapidity region of relativistic 
heavy ion collisions. Assuming that the number baryons of a given kind is given, in first 
approximation, by the Fermi distribution at the freeze out temperature, the relative yields 
of two baryons will be a function of their mass splittings. As a correction to this model, it 
is natural to use the temperature dependent masses in the Fermi distribution (the precise 
sense in which this gives the leading correction will be the subject of a future publication). It 
may be possible then to test our predictions by measuring ratios of the number of different 
baryons species. It is not clear however that this simple freeze out models are accurate 
within 20% so that the relatively small mass shifts can be observed. 

The imaginary part of the pole is shown in Fig. 4. Note that this width can be caused 
not only by the decay of the incoming particle, but by other processes without counterpart 
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at zero temperature, such as the absorbtion of a pion from the thermal background [[IS 
even stable particles can have a finite width at nonzero temperature. For the temperatures 
where we expect xPT to be reliable, say, T < 200 MeV, the width of the nucleons is small 
compared to their masses, and they therefore behave as sharp, well-defined states. 

A hint of deconfinement is found when the residue of the baryon poles is calculated. 
This computation is very similar to the one presented above, the main difference being that 
now that both the octet and the decuplet contribute for all states. For example, the nucleon 
residue changes by an ammount 

27 

^tZ-^ = —H{0) + 15H{-fi), (20) 



where 

^^^^ " 127r2/2 7o ^'^ Uk {x^-u^ + iey 

At T = 200 MeV, for instance, the residue is reduced by 15% in relation to the zero 
temperature value, and decreases very fast with higher temperatures. This is consistent 
with the picture that the baryon poles disappear at high temperature. 



IV. CHIRAL LIMIT 

It is customary to consider the chiral limit where the quark masses vanish when using 
chiral models at finite tempertaure. This limit has the theoretically pleasing feature that 
the only mass scales in the problem are the temperature T and the chiral expansion scale, 
and all quantities have a power series expansion in T. We note here, however, that this limit 
is usefull when T ^ 771,^, which is not of great phenomenological interest given the fact that 
the critical temperature is numerically close to the pion mass. 

Nonetheless, it is worth considering this limit as a consistency check. In the chiral limit 
— ^ 0, our main result Eq. ([T^ ) is 
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This is not a power series in T even in the chiral limit because of the presence of the 
additional mass scale yU (the decuplet-octet mass difference). It becomes a power series in 
T in the limit — >■ 0, in which case the baryon mass shift vanishes. It also becomes a power 
series when — > cxd, corresponding to decoupling the decuplet states; This is the limit which 
corresponds to "pure" HBxPT. In this limit, we obtain 

5TM = T^ + 0(l//i^). (23) 

This is consistent with a general theorem that states that there is no contribution at order 
in the chiral limit. This is also in a agreement with the calculation of Ref. [Q) which used a 
relativistic effective lagrangian for the baryons and concluded that there are no contributions 
to the baryon mass at order in the chiral limit. 
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V. CONCLUSION 



We have computed the temperature dependence of the pole mass of baryons in leading 
order in HBxPT and the expansion for temperatures T < m-Tr. The full dependence 
on the temperature and on the octet-decuplet mass difference is included, as dictated by 
large- iVc and chiral power counting. The shifts in mass were of the order of 20%, at T ~ 150 
MeV, and the changes in width are about the same size. 

The range of temperatures where our results are applicable are limited by both the chiral 
expansion (which breaks down at energies of the order of QCD scales, like the rho mass) and 
by the existence of heavier resonances. The effect of virtual loops of these heavy particles 
at zero temperature are included in the value of the phenomenological coefficients of the 
effective lagrangian. Although such particles are also present in the thermal plasma, they 
are suppressed by large Boltzmann factors at low temperatures, but at higher temperatures, 
they will become important. 

It is important to realize that the large- A^c expansion entered in the present calculation 
in a rather non-essencial way. As long as we decide to include the decuplet as a propagating 
field, consider the octet-decuplet splitting of the order ~ m^, and use a value for the octet- 
decuplet coupling constant consistent with its large-Ac value, our results are essentially 
unaltered. Only a numerically very small term was dropped (since it was 1/Ac suppressed) 
in arriving at ( |l^) . 
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FIGURE CAPTIONS 



Fig. 1- Contribution of order Nml. to the baryon mass. Full lines denote baryons, dashed 
lines with a tick stand for the temperature dependent piece of the meson propagator. 

Fig. 2- Leading non vanishing contribution to the temperature dependent baryon mass. 
The crosses stand for mass insertions. 

Fig. 3- Real part of the mass pole as a function of temperature. 

Fig. 4- Imaginary part of the mass pole as a function of temperature. 
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